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A new variant of the exact Fradkin representation of the Green’s function G c {x, y\gU), defined for 
arbitrary external potential U, is presented. Although this new approach is very similar in spirit to 
that previously derived by Fried and Gabellini, for certain calculations this specific variant, with its 
prescribed approximations, is more readily utilizable. Application of the simplest of these forms is 
made to the A$ 4 theory in four dimensions. As an independent check of these approximate forms, an 
improved version of the Schwinger-DeWitt asymptotic expansion of parametrix function is derived. 

PACS number(s): ll.10.-z, ll.10.Kk, 02.30Mv, 03.65.Db 


I. Introduction 


We here present a new approach for evaluating the (causal) Green’s function G c (x,y\gU) and the related loop¬ 
generating functional L[gU\, for a particle in the presence of an arbitrary external field U{x). For brevity’s sake, the 
simplest case when both the propagating particle and the external field are scalars is considered. The more interesting 
cases, when the particle has spinor structure and/or when it moves in a background gauge or gravitational field, will 
be presented in the next paper of this series. 

The causal Green’s function for a scalar particle in the presence of a scalar background U(x) is a solution of the 
equation])] 

[Px - m - + 9U(x)]G c (x, y\gU) =-S(x-y), (1.1) 

and can be presented in the form 


POO 

G c (x,y\gU)=i dse~ lsm -g(i 

Jo 


,y\gU), 


( 1 . 2 ) 


where g(s\x,y\gU) := e ls ^ +9U ^5 (x — y) is the Parametrix function (with name taken from ref. ||) of the 
problem. Since its relation to both Green’s function and loop-generating functional are very close, we will dedicate 
the major part of this paper to a derivation of different exact and approximate forms of the parametrix function, and 
results from the existing literature will be rephrased in terms of it. There exist two essentially different approaches 
for the construction of the parametrix: the Fradkin forms and the Schwinger-DeWitt asymptotic expansion. The first 
approach will be explained here, while the second is postponed until we reach Section IX. 

The first approach is based on Fradkin’s [3] exact representation 


g(s\x, y\gU ) = e~ 1 ^ dtu ( y+ ti v ) S (x - y - J’ v'j 


(1.3) 


where e l fo dtS v denotes the linkage operator acting on the auxiliary variable v(t) (0 < t < s), and S% := • 

This form, on the one hand, leads directly to a path-integral formulation (Appendix i) , a transition based on 
equivalence relations betw een the functional-linkage and functional integral forms (Appendix |A|) . On the other hand, 
Fradkin‘s expression O) allows a derivation of an entire spectrum of different applications and approximations ( j|, 
||].). The majority of these are eikonal inspired, such as the Phase Averaged Approximation (< Ph >) of | 6 | : 


1 Notation: a := b defines a in terms of b\ dp := Pg '■= Ilf— Metric has signature (2 — D): ry** = 77 ** = 

diag(+l, —1, • ■ •, —1). Also: z := x - y; cr(£) := y + £ 2 ; mi := m 2 — ie\ W M (x) := d^U(: r); H^(x) := d^d^Uix) ; 

We will use the hat symbol for operators that are of continuous nature (differential or integral) ( K ), the hacheck symbol 
for matrices, (7t) and the bold face ) for objects that are both functional operators and matrices. 
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( 1 . 4 ) 


„ <Ph> 


,y\gU) = J dpe ipz+isp2 j dqdPe- i ^~ y ^ p+2 ^e ig i o dtu (q+tP) 


It turns out that, based on the < Ph >-Approximation as an initial step, it is possible to recast Fradkin’s form 
( |Q| ) in terms of Gaussian-weighted integrals over the set of auxiliary 4-vectors {Q„,P„|n = 1,3,5, •••}, with the 
exact representation ( 1- 0): 


g(s\x,y\gU) = J dpe lpz+lsp2 J dqdPe^^^Yl^ J dQ n dP n e~^ p2+Q ^ 

e *afo dtu('q+t(P-2p)-^J2„'i[Pr‘ cos ( 2] f 1 )+Qn sin( a f t )]) 


(1.5) 


where the prime (') denotes summation over odd natural numbers, {Q n ,Pnj are pairs of auxiliary 4-vectors. The 
steps needed to obtain this form are: Taylor-expand the exponential e l9 ^° dtu ( y+ fo w ); use Fourier representations for 
each U- factor and for 5 (x — y — v)\ and evaluate the v- variable linkage operation. This generates exponentials 

with binary products ki ■ kj of the Fourier momenta kj for different [/-factors, and these products are multiplied 
by functions of the difference of corresponding auxiliary (proper) time variables U — tj. One of these functions is 
\ti — tj\; and a suitable separation in terms of its Fourier modes has to be used to allow the factorization of expressions 
corresponding to different U factors. Once that separation is achieved, the reverse (inverse) Fourier transformation 
has to be performed, a nd t he resulting series of factorized terms must be resumed. After some trivial rescalings, the 
resulting expression is (1.5). 

This form provides one with the possibility of approximating the full (exact) expression by taking into account only 
several pairs of variables \Q n ,Pn\- The < Ph > approximation is obtained by neglecting all such {Q ra ,P n } variables 
in the argument U in fll.5|); also, the term — 4sP 2 should be dropped ||. Taking into account only the first pair of 
auxiliary variables {Qi,Pi} in expression (1.5) gives an improved (”< Ph |1 >”) approximation, etc ... It is possible 
to quantify the quality of such approximations, in terms of the quality of corresponding approximations to the | ti — tj \ 
function; and, according to ||, the < Ph\l > may be expected to carry 81% of contribution to the total expression, 
the < Ph\l\3 > 90%, the < Ph\l\3\5 > 93% etc ... 


II. New representation for parametrix 

In the course of searching for an alternate derivation of Q from Fradkin’s form Q (without resorting to a 
perturbative-like expansion and subsequent resummation), we have found a new, exact form, based on the mode 
expansion: 


v(t) p = V p + ^2,\Qn cos (nut) + P p sin(nwf)] 


( 2 . 1 ) 


71=1 


where w = —. The main criteria in developing such an expansion is (the ad hoc) requirement that the delta-function 
restriction dtv{t) = z constrains only the amplitude of the zero-frequency mode V (forcing it to be = j), thus 
leaving amplitudes of the oscillating modes unconstrained. The properly normalized basis on the interval [0, s] is {( f > n } 

= jco(t) = yi; Cn(t) = cos(nwt); s n (t) = sin(naA); (n > 1) j, where the scalar product is just the integral 
over that interval. 

With such a basis the resolution of the identity on the same interval is provided: 


5{t (j>x{t)*(j)x{t') = - + - ^ cos (nu;(t - t'))\ 


( 2 . 2 ) 


71=1 


and, requiring that = S(t — t')5 p , the corresponding mode-decomposition of the functional derivative^ : 


2 For the general decomposition v{t) p = V g <f>{t)\ where {0 (/)a} is an ortho-normal basis; the corresponding functional 

derivative is given by m Ea W£- 
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(2.3) 


Sv(ty 


1 d 
s dVv 


n= 1 


- ^ < cos (nut) 


dQt 


sin(nujt) 


OPn 


Using (2.3), one can calculate some necessary expressions: 

S 2 


dt 


5v{ty 


1 d 2 
~sW I 


r) ^ 

-v 


d 2 

dQl 


d 2 1 

dPU 


ft S °° 1 

/ dt'v{t') = Vt ^-— {Q n sin {nut) + P n [ 1 — cos(nwt)]} 

Jo 2tt ^ 7Z 

f dtv(t) = Vs ; 

Jo 


(2.4) 


and evaluate the U-linkage operation (using eq. ( |A2| )), to obtain the new form for the parametrix CD : 


g(s|.T,J/|ffU) = 


(47ris) J 


n 


( agj efj J 


7 Jq d.tU(y+z- 1+ 2 ^ 5Z ^=i sin(nwt)+P„[l—cos(nijjt)]}) 


{Qn},{P„}-0 


(2.5) 


This form is an exact representation of the Fradkin parametrix. Note that the characteristic ” freq uencies” for 
modes in auxiliary space [0, s] employed here are even multiples of ((, as opposite to the FG form (L5), where the 
characteristic frequencies are odd multiples of that value. Also, here the variables V, {P m ,Qm}’s are just Fourier 
mode of the auxiliary variable v(t), while in form (|l^) no such simple connection between v(t) and their { P m , Q m }’s 
is known. 

In the following Section we will develop several approximate forms to (2.5), as well as one suitable reformulation of 
that exact result. That reformulation, in it’s turn, then will give rise to the some other approximate forms. 


III. < 0 > Approximation 


We are now able to develop the hierarchy of approximate forms (in the same spirit as those of the < P/i|l|3| ■ • • >). 
The basic one is obtained by neglecting all oscillatory modes (Q n ‘s and P n l s) in the exponential argument of U in eq. 

g <0> (s\x , y\gU) = —L-^e-^Jo dtu Hi)) ( 3 . 1 ) 

(47ns) 2 

where a (a) = y + az. We will call this the < 0 >- Approximation. It has a simple meaning in coordinate 
space: propagation is replaced by motion along the straight line between points x and y. In other words, the < 0 >- 
approximation ignores the ’’structure” of the potential U(x ), paying no attention to it’s ’’valleys” and ’’hills”. In terms 
of a propagating particle, this is the extreme high-energy (UV) approximation; in terms of a background potential, it 
is requirement that all its spatial derivatives are ’’small” (i.e. \d^ ■ ■ -d^JU\ « where the characteristic length 
L is ~ s < v 2 > 1//2 ). 

It is worth noticing that the corresponding Green’s function has the form of a free particle’s Green’s function with 
a path-dependent mass: 


G<" > (x,y\gU) = - 


ds- 


(47t is)' 


a -isM(x,y) 2 _ - 


= G c {x,y\Q) 


m 2 — (x,y) 2 


(3.2) 


where M(x,y) 2 = m 2 — g f Q daU(a(a)). From this form are evident the properties of symmetry, under x <-> y; and 
of the correct free-particle limit, U —> 0. 

The < 0 >-approximation is not same as the < Ph >. This can be seen by comparing a transformed (as closely as 
possible to the < Ph >-form; see Appendix |6) expression for < 0 >: 
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g < 0 > (x,y\gU) = j dpe ipz+isp2 J dqdPe- iq{p+ 2 p) e lg ^ dtu (y+^+tp) ^ 


( 3 . 3 ) 


with: 


g <ph> {x,y\gU) = J dpe ipz+isp2 J dqdPe~ iq{p+2p) e ig ^ dtu( y +q+tp) 


(3.4) 


The two forms are very similar; the only difference between them is in the argument of the potential U: instead of 
[y — q — tP] (in < Ph >) we have [y — f <7 — tP] (in < 0 >). 


While the < 0 >-approximation has the free-particle-like coordinate form (3.2) for the Green’s function, no such sim¬ 
plification can be made for the < Ph >-approximation. But, in momentum space picture, the < Ph >-approximation 
has a nice interpretation while the < 0 > remains nontransparent. Thus one may consider the < 0 >-approximation 
better suited for coordinate-space applications, while for diagrammatic applications, the < Ph > seems better adapted. 

One more property (besides x <-> y symmetry and the correct free-particle limit) that is shared by both the < 0 >- 
and < Ph >-approximation is their exactness for the plane-wave (laser-like) background U(x) = <j>( k.x) (with fc 2 =0). 
The final form of the exact Green’s function of this field is the given by the same expression, Q3.2| ), which is only the 
< 0 >-approximation for other (non-laser) backgrounds. Also, it is worth noticing that in the family of approximations 
based on forms ( |3.3| ) and (p.4[), where the argument of the potential U is replaced by [y — qF\{a) — sPF 2 (a)], the 
laser-background case again simplifies to the explicit < 0 >-form (3.2). 


IV. The loop generating functional L[gU] 

Another important object (besides the Green’s functions) for functional field-theoretic calculations is the loop 
generating functional: 


L [gU] ■= —-Trdet 


(P 2 -m 2 _+ gU x ) 
(■ P 2 -m 2 _) 


= iTrdet(G c [yC/]G c [0]^ 1 ) 


(4.1) 


It can be expressed in the integral form L\gU] = \ f () dXTr (G c [\gU]gU) = | f Q dX f dxG c (x,x\\gU)gU(x). 
In the < 0 >-Approximation ('[i.2|) the Green’s function has a diagonal component of amount 


G<°>{x,x\gU) = - 


ds- 


(47 vis)' 


0 -is(m 2 _-gU) _ 


{A-k)- 


and one finds 


■0 m 2 _-gU )" 1 I' ( 1 - ^ 


(4.2) 


-<o> 


[gu} = 


d\ / dx- 


2(477)"? 


(47r) J 


-{rn 2 - XgU^-'-gUT 1 - 


D 


dx 


(m 2 — gU ) 2 — 


(4.3) 


This is a divergent expression for even P, with T (— -2) expressing the UV-divergences. In what follows, we will 
assume that T ( — 4^) is regularized in some implicit wayjj, allowing explicit evaluation of < 0 >-approximate forms of 


some simple field theoretical models. 


Note that for D = 2 and D = 4, expression (4.3) simplifies further to 


3 One such regularization follows from the definition of the F-function, by replacing it’s UV-limit (p = 0) with the finite value 
A/ 2 A" 2 : 


r - 


-§H 


M 2 A - 2 


dpp 2 x e p 


j 

J A 


M 2 A ” 2 


dpp 


where M is some fixed mass scale. Therefore, r (— ^) ~ A 1 


2 

D 



(4.4) 
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L <0> [gU} D=2 = J dxU(x) 

L <0> [gU ] D -4 = J dx[-2m 2 gU(x) + g 2 U[x) 2 } (4.5) 


V. Application to the A<p 4 -model 


To illustrate the < 0 >-approximate form, consider the Xp 4 model in Minkowski 1 + d-dimensional space. Start 
from an action 


S= dx 


\[dp) 2 - 


(5.1) 


which is (after evaluating the EOM of the field ip) equivalent to the standard action with a potential term V (p) = 
— i m 2 if 2 — jf 4 , where A = 2g 2 /p 2 . This particular formulation is chosen in preparation for subsequent steps where 
it is important to have the action with all terms at most quadratic with respect to the any specific field (for example, 
the interaction term gp 2 ip is quadratic with respect to the ip field and linear with respect to the ip field). The field ip 
may be thought of as the classical form of a composite ip 2 . 

The generating functional Z[J,j) := J{D(pDip]exp { — iS[ip, ip] — i J dx[pJ + ipj]} can be transformed, following 
steps outlined in Q, into the form^j 


Z[J. j ] = g h f f i 


(5.2) 


where G c [gU] := [d 2 + m 2 _ — gU) 1 is the causal Green’s function for the field p in an external potential [/, and 
L\gU] := ilndet (^G c [gU]^ is the corresponding vacuum-loop generating functional. 

In this section, the the 2n-point quantum correlation functions 


r 2 n(aq, • • -X2n) ■= < [p[xi) ■ ■ ■ ip(x 2 n))+ >= 


1 


■2 n 


5 2n Z 


Z[ 0,0] Sj[xi) ■ ■ ■ SJ[x 2 n) 


J,j =o 


= i n [G c (xi, x 2 \2gi8j) ■■■ G c [x 2n -i, x 2n \2giSj) + permutations] Z[ 


> 3=0 


(5.3) 


will be evaluated in the < 0 >-approximation. Here Z\j\ ■= Z[0, j]/Z[0, 0] = Z[0,0] 4 exp [L[2gi5j]} exp j ^ f j 2 |. 
Only the first of these terms will be observed, having all others related to this one by a permutations of coordinate 
arguments {xi, ■ • ■, x 2 n } ■ 

For D = 4 use H) and dA4| ) to obtain 




(5.4) 


where a := r R 2 J X . From here, use (El) 


r 2 n [x\ , 


• • • X2, 


0 = H)"II 


ds 


4j- (X2j-\~ x 2j) 2 -im 2 _Y, i Sj 


3=1 1 


3 ( 47 risj) % 


r<0> 


-2 gS^Sj / daS[y - x 2j - a(x 2 j-i - x 2j )) 
j J ° 

The Z <0> factor can be represented as a product of two factors K\K 2 where 


(5.5) 


4 An alternate form is Z[J, j] = e ^ 1 d l e L M~P I J-a c [g^]J 


1p = ~- 
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(5.6) 


Jd ■ = g( i^rj) E"=l E" = l SlSr JgdaJgd/3S(x2l—X 2r +a(x2l-l—X2l)-l3(x2r-l-X2 r )) 

I<2 := e -*T^5:”=i S| 


The exponential in A'i has two sources of singularities: First, in ’’diagonal” terms l = r, the strongly UV-singular 
integrals J* da df3S£>((a — /3)(x 2 r-i~X 2 r)) is present. It can be transformed to the form f dk ^(^ 2 r-i ) _1 sin^ 2 r-i^ , 

where £ 2r _i — \k\x 2 r-i~ ^ 2 r), that it is well-behaved only in one direction (parallel to the four-vector (x 2 r-i—X 2 r ) fJ- )- 
This gives an overall UV-behaviour ~ A 3 , where A is UV-cut-off parameter, the same as in (4.4). 

A second source of singularities in K\ are nondiagonal terms in case of intersecting linear trajectories; such 
trajectories occur whenever there exist a and (3 such that X 21 + a(x 2 i-i — X 21 ) = X 2 r + P(x 2 r -i — % 2 r )■ 
Their contribution da J d^Sn(x 2 i + a(x 2 i-i — X 21 ) — X 2 r — (3(x2 r -i — % 2 r)) can be transformed into the form 
/ dk (£ 2 r— 1 ^ 2 /— 1 ) _ 1 sin^ 2 r-i sin £ 21 -1 exp {| ( X 21 + X 21-1 — X 2 r — X 2 r-i ) • k}. From that expression one can infer that 
two fc-space directions are convergent, and in one more the integration is made convergent by a fast-oscillating 
exponential factor. The remaining unregularized direction gives an overall A 1 UV-behaviour for such intersecting 
trajectories. 

Then, since a ~ T(—2) ~ A 4 , terms in the exponential of K\ will display one of three UV-behaviours: A 3 ~ 4 = A -1 
(diagonal singular terms), A 1-4 = A ~ 3 (nondiagonal singular terms) or A ~ 4 (nondiagonal nonsingular terms, i.e. 
nonintersecting trajectories). In this way, in the A —> 00 limit, one obtains K\ = 1. 

Notice that in the A —> 00 limit, the value of bare-coupling constant A does not play any role, since it cancels: 
A = —> — r ^ 2 ) ~ A -4 —> 0. Only way to make A nonzero (and finite) in the A —> 00 limit is to accept peculiar 

behaviour of bare coupling A = — 2) ~ A ~ 4 —> 0, which is comletely opposite to the usual behaviour of bare 


1 - 


*r(- 2 ) 

8,2 


couplings in the Abelian models. 

Concerning the factor A' 2) there are no other (than a) singularly behaving quantities in the exponential, and one 
can take the A —» 00 limit right away. The result is A ' 2 = e + lm2 E"=i s ‘. Or, since is not necessary to take that limit 
at this point, one can join it to the bare-mass terms in the expression for the r 2n , having, in effect, replaced the 
bare-mass m 2 with M 2 := m 2 + = jh^m 2 in all subsequent expressions. 

Then the 2n-point correlator factorizes: 


^2n (x\ , X2n) + 11 G c {X2j — 1 , X2j 10 ) 

J =1 


m 2 — >M 2 


(5.7) 


into product of a free-particle Green functions. 

Alternatively, one can evaluate the generating functional in naive T(—2) —> oo limit to get 


?<o> 


[J, 0 ] Z <u> [ 0 , 0 ]e 


4fJG< 0> [-aM 2 ]j 


= Z 


<o> 


[ 0 , 0 ] e 


U JGc 




(5.8) 


i.e. such a limit produces a free theo ry of fields of mass M. Notice that it is not necessary to take the < 0 >- 
approximate form for G c [gU] in (|5.8| ): simplification in the K\ factor in Z <0> [j] suffices to validate the above 
expression with the exact G c [gU]’s. 

The disappearance of the effective interaction in the UV-limit of tp 4 theory in four dimensions is the well-known 
phenomena of triviality, established or indicated in several different frameworks (Euclidean constructive field theory: 
[ |l2|| - Jl5|; L attice calculations and ’’inspired” perturbative-renorm-group extrapolations: 0- n ; Functional ap- 
proach: |l(i[). Here, we have presented a simple analytic derivation of triviality. However, this method is not uniformly 
suitable for all number of dimensions: for example, in D = 2 the quadratic linkage operation is absent, leading to the 
absence of UV-singular quantities in denominators (that were o f su ch importance in the D = 4 case), and thus to the 
non-triviality of the theory. In D ^ 2,4, the linkage functional (T4) is either (in case of even D) a polynomial of order 
higher than two, or a nonpolynomial (irrational) function (in the case of odd D) of gU ; both circumstances generate 
linkage operators e L ^ 2l9Sj ^ that are neither shift-operators nor Gaussian linkages, preventing further evaluation in the 
absence of additional approximate steps. 

The more detailed discussion, concerning more uniform regularization of sources UV-singularities, how this influence 
the derived trivial behaviour, and corrections due to higher approximations (using < L >,< Q > and approximations 
obtained using the Schwinger-DeWitt’s approach; all these are introduced in next few sections) will be presented 
elsewhere 0 . 
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VI. Beyond the < 0 >-Approximation 


How may one proceed b eyon d the < 0 >-approximation? The goal here is to obtain a sequentially better ap¬ 
proximation to expression (2.5). Imitating the < P/i|l|3|--- >-hierarchy does not help much, since the potential 
U(x ) is in general a nonlinear function of x : for example, retaining only n = 1 level variables Q\ and P\ (the 
< 0|1 >-Approximation) gives: 


„<!>/ 


_lif a _i_ a \ . . 

(s\x,y\gU) = g(s\x,y\0)e S l / e i 9 s Jo daU(<r(a)+^:{Qi sin(27wra) + Pi[l-cos(2rMra:)]}) 


( 6 . 1 ) 


This expression can be evaluated analytically only in some special cases: e.g. if the potential U(x) is a quadratic 
polynomial of its arguments, or if it represents a laser field (when it reduces to < 0 >-approximation). For all other 
cases one has to resort to some kind of further approximation, such as a Gaussian approximation. The same situation 
holds for < 0| 112 > and all higher approximations. 

This indicates that it is perhaps not beneficial to insist on an application of the hierarchy < 011121 • • • >; but, 
instead, to proceed by an alternate route. Let us, in < 0 >, ’’remove” the variable V from v(t), while keeping all 
other components of v(t) grouped into a new variable u[t) v(t) — V. The variational derivative S v then splits 
into 5 v (ty = \dva- + &' u t t y Here the is a constrained variational derivative that takes into account the 

constraint ff dtu(t ) M = 0; its action on u{t') v is t'), where n(f, t ') = [5(t — t') — i] is a projector on 

the u-subspace. In this way, one can write 5' u = n5 u in terms of the unconstrained variational derivative S u . 

Then, the exact expression for the parametrix can be reformulated (Appendix [|]) as 


9 (s\x,y\gU) 


(47ns) : 


-iflg-Go 3 S 'l J9J0 dtU ( a (j)+Jo u ) 


u —>0 


( 6 . 2 ) 


While the V-integration (linkage operation) was easily performed in an explicit manner (using identity ( |A2| )), the 
u-integration cannot be done in an exact way. Instead, one can assume that the u variable is a ’’small” quantity^, 
thus generating a Taylor expansion of the U [cr(t/s) + in the u variable, and keeping only some of the lowest 

order terms. In what follows, derivation of the two such lowest order approximations is defined by retaining only up 
to 0{v}) (giving the linear in u, < L >-Approximation) and 0(u 2 ) (quadratic in it, the < Q >- Approximation) 
terms. 


VII. < L >-Approximation 


The first of these, the < L >-approximate form for the parametrix is a product of the < 0 >-approximation 
parametrix and a linkage operation factor 


= gGo “ 2 -4(/o °) 2 


u—> 0 


e -GoXVJo 

where a^t) := gs ft daW^ (<r(a)) and W tJi := d^U. Using formulas from Appendix £j], one can obtain: 

/»S /*1 nCX. 

/ a 2 = 2s 3 g 2 / da d/3/3W ll (a(a))W IJ ,(a(P)) 

Jo Jo Jo 


f a) = s 3 g 2 f daa f d^W^a^W^ 

Jo J Jo Jo 


(7.1) 


(7.2) 


(7.3) 


so that the final expression for the < L >-approximation is: 


5 A similar approach is called the string inspired formalism by Bern-Kosower [jioj. |ll[: the zero-mode variable is separated 
and the expansion over the oscillating modes is performed, allowing the efficient resummation of many parts of perturbative 
Feynman graphs 
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g <L> (s|s, y\gU) = g <0> (s|®, y\gU)e 2is3g2 So **(i-a)w„("(°))/“ dppw^m 
At z = 0 it simplifies to: 


9 > (s\x,x\gU) = 


(4-7TJS)' 


JsgU(x)+i^-W 2 /x) 


(7.4) 


(7.5) 


VIII. Q >-Approximation 

In the case of the < Q >-Approximation, one retain up to quadratic terms (in the u variable) in the Taylor expansion 
of U. The parametrix is again a product of the < 0 > form and a factor with the functional operations: 


g */o g * Jo u ' a ~^~ 2 Jo Jo 


_ - iTr ln(l-2II-b)+i/ / a-(l-2n-b) _1 -n-a 


u —>0 


( 8 . 1 ) 


where bp V (t,t') := gs fmaxt.t' da/H^v (cr(a)), H{x) := (9® dU)(x), and eq. (|A4j) has been used. 

The determinantal term tV[gJ7] := — iTrln(l — 211 • b) = 2 n Tr[(II • b) ra ] can be evaluated term by term, 

but no general closed form can be obtained. The first two terms in that sum are (Appendix ^): 


N[gU ]i = +gs 2 f daa( 1 — a)(d 2 U)(a(a)) 

Jo 

pi pd. 

N[gU] 2 = 2g 2 s A / da(l - a) 2 H^(a(a)) d/3f3 2 H^(a((3)) 

Jo Jo 


( 8 . 2 ) 


The value of N[gU] at z = 0 can be obtained explicitly, using a different, more direct approach (Appendix |g|), with 
the result: 


( S y^)\ 


= t (Mwu.) = 


1 / sin i 

N\gU}\ _ n =-Trln _ . — , 

2 ( s^2pH ) 4 ^ 771 ( 2771 )! 

gs 2 . g 2 s 4 4g 3 s 6 g 4 s 8 16g 5 s 10 

= V £lW + Ho C2{x) + 28 K £a(l) + 47ffi £l(l) + 4CT7re £s(l) + 


(8.3) 


where C m (x) := Tri(7i m ) and Tr^ is the trace over Lorentz indices^. 

In a similar way, for z / 0, one can obtain the Tith-order terms in Taylor expansion of the expression \ — 
if f a ■ (1 — 211 • b ) _1 • II a. The first two are (Appendix ||): 

Xo ■= i J J a, ■ II • a = i2s 3 g 2 J da(l — a)W ll {a{a)) J d/3pW‘ 1 (a(/3)) 

Xi := 2 i J f a ■ Tl ■ h ■ II ■ a 

= i2s 5 g 3 ( da f d/3 f d 7 (min(o, 7 ) - cry)(min(/3, 7 ) - f3x)W, J .{(j{a))'H' 1 ''{a{i))W„{a//3)) (8.4) 

J 0 Jo Jo 

Notice that retaining only \o is giving the < L >-approximation. 

It is possible to obtain an closed exact form for \ at z = 0 (Appendix |): 


Xl 2=0 = 


isg 


ban f s\! 


W ■ 


n 


\ 


- 1 


w 


(8.5) 


For example, several first members are: 


example: TriR := R a 


















( 8 . 6 ) 


XoUo = i S ~^W 2 (x) 

XI Uo =i^-(W-n-W)(x) 


X 2 \ z =o = l 


60 

17 

5040 


g 4 s 7 (W-H A ■ W) 


The closed exact form of < Q >-parametrix is then: 


_ \ f | ^tan (s\J ^ j ^ 

g<-> { s\x,x\ g U)=i exp ^-sgW ■ ± ^ 1 ■ 


IT 


(4-7Tis)' 


3is9t / + 9^£ 1+i £J_ W 2 + £ « £2+i£ J_ (lv . H . w)+ ^= £3+i 17 9 4 s 7 (ly . K 2 . w)+ |* £4(;i;) + .. 


(8.7) 


IX. The Schwinger-DeWitt Asymptotic Expansion of the Parametrix 

The Schwinger-DeWitt Asymptotic Expansion || method allows an arbitrarily precise approximation of the 
parametrix g(x, y\gU). The method is known under many different names (heat-kernel method, high-temperature 
asymptotic expansion, etc.) and was applyed in various problems (spectral geometry of manifolds with boundaries 
||, propagation of fields in general relativity |fo|- j2fj], etc.). We are presenting here its simplest variant with the 
slight improvement allowed buy one-component nature of propagating held. 

Starting from the Initial Value Problem (IVP) for the parametrix g(x,y\gU): 

- id s g(s\x, y\gU) = [-d 2 + gU(x)]g(s\x, y\gU) 

g(s = Ojac, y\gU) = 5{x - y) (9.1) 

one may solve it in two stages: Firstly for the free case U = 0, with the solution go(s|x, y) := g(s\x, y|0) = ——n-e _I ^. 

(47rs)"3" 

Secondly, for nonzero potential U, introduce the ansatz g = go ■ h. The IVP for the reduced Parametrix h is then: 

- i (d s + ■ d x ^j h(s\x,y\gU) = [-d 2 + gU{x))h{s\x,y\gU) 

h(s = 0\x,y\gU) = 1 (9.2) 

Although it is not possible to solve the first of these equations exactly, the method allows one to reach (in principle) 
any desired precision. Its idea is based on the notion that the differential operators on the LHS and RHS of IVP ([h|) 
are homogeneous of degrees —1 and 0 respectively. Then, for the formal Taylor series h(s) = 1 + i(i s ) n h n (this 
is an asymptotic series), one will find a well-separated system (hierarchy) of recursive equations for the functions h n : 

(z ■ d x + l)hi = gU 
(z ■ d x + 2 )h 2 = ~d 2 hi 
(z ■ d x + 3 )h 3 = -d 2 ,h 2 + gUh.2 
(z ■ d x + 4 ) 6,4 = ~d 2 h 3 + gUh 3 
(z ■ d x + 5)65 = -d 2 h,4 + gUhi 

(9.3) 

These can be solved in sequence, beginning with hi, and continued to any finite order. Regularities will ap¬ 
pear, suggesting that a better (in terms of computational costs) starting point is to use the exponential form 
h = exp{ VOL-, (is) n k n |. One can always switch from one form to another, using algebraic relations which con¬ 
nect the two hierarchies of the coefficient functions {h n (x,y)} and { k n (x , 9)0 The hierarchy of equations for the 
coefficient functions k n is: 


hi — fci, 62 — 62 + 2 , 63 — 63 + kik 2 T 64 — ki -t- kik 3 -t- ^ k 2 -t- 2 61 k 2 -t- 24 61 >* * • 
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(z-d x + l)fci = gU 

(z-d x + 2 )k 2 = -dlki 

{z-d x + 3 )k 3 = -azk 2 - (dzfci) 2 

(2 ■ + 4)fc 4 = -<9 2 fc 3 - 2d x k 1 d x k 2 

(z ■ d x + 5)fc 5 = - ( d x k 2 ) 2 - 2d x k 1 d x k 3 

(9.4) 

One can resolve this hierarchy iteratively up to any desired level of accuracy 0(s n ). The first several terms are 
(Appendix |k|): 


ki(x,y) =g daU(a(a)) 

Jo 

k 2 (x,y) = —g I daa(l — a)(d 2 U)(a(a)) 

J 0 

1 f 1 

h(x,y)=g- J daa 2 (l — a) 2 (d 4 U)(a(a))+ 

pi pCt 

-2g 2 / da{ 1 - a)W(a(a)) / d00W((r(0)) 

J 0 J 0 

k 4 {x,y) = ~g J daa 3 (l — a) 2 (d 6 U)(a(a))+ 

pi pOi 

+g 2 daa(l — a) 2 (d 3 U)(a(a)) / d00W(a(0))+ 
Jo Jo 

pi pCX. 

+2 g 2 da(l-a) 2 H^(a(a)) d00 2 H^ (<r(/?))+ 

Jo Jo 

+g 2 [ — (l~a) 2 W(a(a)) [ d/3/3 3 (d 3 U)(a(0))+ 
Jo a Jo 


+g 2 [ — (1 - a)W(a(a)) f d/3/3 2 [2a - 0(1 + a)\(d 3 U)(a(0))+ 

J 0 a Jo 

+g 2 f daaW(a(a)) f d00(l - 0) 2 (d 3 U)(a(0)) + ■ ■ ■ 

Jo J a 


Although straight-forward, the calculations become very tedious in higher orders. 

The (exponential-) asymptotic expansion of the parametrix has the approximate form: 


(9.5) 


g(s\x,y\gU) 


^ c -i£ c isQ fj daU(cr(a))+s 2 gfg daa(l-a)(d 2 U)(cr(a)) 

(4m s)^ 

e ~is 3 [^gfg daa 2 (l-a) 2 (d 4 U)(a(a))-2g 2 Jq da(l-ct)W(cr(a)) /“ d/3/3W(<T(/3))]+s 4 fe 4 + '" 


(9.6) 


and one should note the agreement of the U -dependent exponential in that expansion ( |9.6| ) with terms appearing in 
previous approximate forms: for example, the sg term comes from the < 0 >-approximation; the s 3 g 2 term comes from 
the < L >-approximation; the s 2 g and s 4 <? 2 terms come from the < Q >-approximation (N 2 [gU] 4 and N 2 [gU] 2 ), as 
will all higher terms that depend on second derivatives matrix TL := (d<S>dU) only (they will be of order ( s 2 g) n ~ s 2 ™); 
the < Q > correction terms for < L >, of the type W ■ H n ■ W (n > 0) are not visible here, since their order is 
s 3 g 2 (s 2 g) n ~ s 3+2n . The term s 3 g is proportional to the fourth derivative (hyper)matrix (d ® 4 [/), and cannot come 
from either <L>or<Q> approximations. Same is with (9® 6 [/) (~ s 4 g) and higher derivative terms. 

The diagonal element (in coordinate space) of the parametrix (i.e. z = 0, <r(a) = x) takes form: 


g(s\x,x\gU) 


1 is g U(x)+^L(d 2 U)(x)-i4Sd i U(x)+iA?J-W(x) 2 

D ^ 

(4ms)~ 

e -^{d 6 U)(x) + ^§J(d 3 U){x)W(x)+^H^{x)n^{x)+- 


(9.7) 


This form can be useful for the evaluation of the functional determinant that represents vacuum loops in quantum 
field theory. 

As a last point in this section, note that is possible to obtain the full expression for the collection Fi(z|s) of the 
first order in g terms in the exponential k = ln/i (i.e. h = e9 F i( z ls)+0(s )). The result is (Appendix |]) 
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(9.8) 


pi pi °° 

Fi(z\s) = is / dae~ lsa ^ 1 ^ a " ,a y U(y + az) = is / da ^ [— isa(l — ce)dy] n U(y + az) 

Jo JO _ _ n 

It’s diagonal (z = 0) value is 


n —0 


^i( 0 |s) = is 


n —0 


(2n + l)H V 2 


--dy) U(y)=isU+-d>U-—d*U-—d°U + 


60 


840 


which is in full agreement with previous derivations. 


(9.9) 


X. Conclusion 

In this paper we have presented a new form of Fradkin’s parametrix, as well as a set of reasonable approximations. 

The first of them, the < 0 >-approximation, is eikonal in spirit, in the sense that it takes into account only the 
extreme high-energy behaviour of the propagating particle. It is exact in the case of a laser-like external field. That 
property will be shown (in the next paper) to be true in more complex cases (spinor particle, gauge and gravitational 
external fields). Since the < 0 >-form has a simple coordinate representation, we did employ it in calculation of full 
(nonperturbative) quantum correlation functions for A<p 4 1 zj =4 theory. The effective interaction vanishes, displaying 
the ’’triviality” of theory. 

Beyond the < 0 >-approximation, two directions were explored: first, < L >- and < Q >- approximation that are 
first two corrections to < 0 >, derived in the spirit of a semiclassical approach; a difference from the semiclassical 
approach is that fluctuations are not taken around the classical trajectories of the particle, but around eikonal ones 
(straight-linear trajectories). Another approach is the Schwinger-DeWitt’s asymptotic expansion method (heat kernel 
method), which we have adapted and made explicit. The validity of its results is established by a comparison with 
< Q\L\Q > results. The majority of longer derivations have been placed into Appendices. 

The author is indebted to H.M.Fried for generous help during both the research and in preparation of this article. 
Thanks are also due to G.Guralnik, K.Kang, A.Jevicki and P. Neskovic for useful conversations. 

This work was suported in part by the U.S. Department of Energy under Contract No. DE-FG02-91ER40688-Task 
A. 


APPENDIX A: Some useful linkage and functional identities 


n —ictd(, 


m 


v = dV ’ „ e~%) 
v^o J ( 47 na) 2 


(Al) 


e~ isd v5 D {V - a)f(V) 


v^o 


(47ts) ; 


-e f{a) 


(A2) 


e -i/o 4*%f(v) = [ [Dv]N(s)e-iJo dtv ^ 2 F[v] 

v^O J 


(A3) 


where N(s) := f Dr)e l fo v = (pDve *Jo dtv (t) ^ 


e J 5 u A-5u e J a-u+f u-B- 


u—> 0 


= det(l - 4AB)~hJ a '^™ A c 


(A4) 
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APPENDIX B: Functional form of Green’s function 


It is possible (using the identity ([A3|) ) to write ( 0 ) as a functional integral over the auxiliary variable v(t ): 

G(x,y\gU) c =iJ dse~ ism - J [Vv] N(s)e~i^ dtv(t) 2 e if° dt g u(y+f^v) 6 ^ _ y _ J \ ( B1 ) 

where N(s) is a normalization constant fixed by t he re lation N(s)~ 1 = f [Dv] e~i Jo dtv W . 

One can transform the functional-integral form (Bl) further into a more-heuristically transparent form, 


G(x,y\gU) c = i dse 


*( pi - p2 -^W)^(l( 0 ) _ y)6(X(s) - a;) 


(B2) 


by introducing an auxiliary variable X(t). Impose on it the constraint X(t) = y + fj £ v, which is equivalent to the 
IVP X(t) = v(t),X(0) = y. The idea is to solve that IVP for v (i.e. v(t) = X(t)). Then one can perform a change of 
variables: v —> X by including the identity (see Appendix Q) 


1 = / DXS 


X(t)-y - / v 
J o 


= / DXS 


X{t)-v(t) S(X(0)-y) 


(Bl), 


poo n 

G(x, y\gU) c = i J dse- ism - J [DX] N{s)S{X( 0) - y)S(X{s - e) - x)e^o~ e dt g u(x(t)) s 


(B3) 


(B4) 


where the expression S is: 


S= [Dv\e-*f° dtv{t)2 6\x(t)-v(t)] = [DPDv]e~if° dtv(t) 2 e ~if° dt[PX-Pv] = h Dp ^ e -iJ’dtPX+iJ‘ 


dtP 


(B5) 


Finally, one may assemble all terms to get (p32|). 

Its heuristic interpretation is very simple; that of a path-integral amplitude for a particle with Hamiltonian 
H(X,P) = P 2 + U(X), which propagates between initial position y and final position x. Since the particle is 
relativistic, the reparametrisational symmetry of its world line must be fixed, i.e. one should integrate over all possi¬ 
ble choices of world-line length s (between these two points), with the ’’weight function” e~ lsm determined by the 
mass of the particle. 


APPENDIX C: Change of variables 


The proof of the identity: 

f DXS 


X(t) — y — 


= / DXS 


X(t)-v(t) S(X(0)-y) 


(Cl) 


may be given as follows. Start with the left hand side (LHS) of ( |cT| ) , discretize it to 
LHS := lim [ dX^dX\ ■ dX^StX o — y)S(X i — y — Atvo) 

N—>oo J 


S(X 2 - y - A t[v 0 + Ui]) • • • S(X N - y - A t[v 0 -\ -1- vjv_i]) 


(C2) 


where At = e =:= and X m := X (mAt) and v m := v (mAt). Note that the very first (5-function does not have a 
u-variable in its argument. Then replace y —> Xq in the second (5-function (using first (5-function), replace vq —> Xl J t X ° 
in the third (5-function (now using the second (5-function), etc ... As a result, one obtains 


lim (At)^- 1 / dX 0 dX! ■ dX N S(X 0 - y)S ( — 
N^oo J V 


Xi - X, 


o 


A 2 - Ad 
St 


-V! •••$ 


St 

Xn — Xn-i 


St 


- V 0 

VN-l 


and then remove the regularization (N —> oo). The final result is 


DXS 


X(t)-v(t) <5(A'(0) — y) 


(C3) 

(C4) 


where the limjv-Kxj(A t) N 1 is included in the new measure DX. 
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APPENDIX D: < 0 >-form versus < Ph >-form 


To transform the < 0 >-form <0 to ( |3.3| ), introduce the auxiliary variable P := -z so that 


Q (s\x,y\gU) = 


(47ris)“ 


fo = 


= i{-ins)% j dPS(z - sP)e-^- +i9 -fo dtu ^+ tp ) = 

= i(_ i7rs )T J d p dpe ip(^P)-^f-+i9f 0 a dtU(y+tP) 

Then represent i(—ins)^ as f dqe* 9 , 

g <0> (s\x,y\gU) = J dPdpdqe ipz - ispP -^ +L ° q2+ia ^ dtu(y+tP) = 

= J dPdpdqe ipz+isp2 -* s( ' F+2p ) 2+ ° q2+i9 So dtu (v+tP)^ 

and shift g—>g+is(.P + 2 p) to obtain 

g <0> (s\x,y\gU) = J dPd P dqe ipz+isp2+ ^ q2+iq{p+2p)+i9 1° dtu{y+tp) = 
— J dPdpdq e ipz+isp2+iq ( p+2p+ i q ) +i9 So d * u (y+t p ) _ 

_ J dPdpdqe ipz+isp2 ~ iq ( p + 2p )+ i 9 fo dtU(y+tP+±q) 

where the last two ste ps w ere shifts P —> P — -q and q —> —q. 

That is expression (T3), suitable for comparation with < Ph >-form 


(Dl) 


(D2) 


(D3) 


APPENDIX E: Proof of projected form 


The easiest way to prove the form (622) is to notice that constrained variable II • u = u — uq (uq := ^ Jf u) is given 


by the expansion (§T|) as n • u{t) = E)T=i [Q£ cos(nwf) + P% sin(nurt)]. Then ff S’l = f ESEi and 


f 0 u = tu 0 + ^ E“=i \ {Qn sin (nut) + P n [l - cos(nwf)]}, giving: 


g(s\x,y\gU) = 


n 


if a 2 , a 2 1 
\ a Qii T ap A J 


(47ris) = n=1 

pa Jo dtu ( < 7 (i) +tu ° + ^ E“=1 7;{Qn sin(no;t)+P n [1—cos(nu;i)]}) 


{Qn},{Pn},U 0 ^0 


(El) 


Sinc e there is no linkage operation with respect to uq, one can directly proceed to the uq —>► 0 limit, reproducing the 
(2.5) form. 


APPENDIX F: Evaluation of N [gUf , 


a — max 


Start from ( |J10| ): 

(II ■ b) pt UT = gs f doth 9 v {a{a)) 

Jo 

where H p v := (d p d v U). Then 

N [gUf = Tr(II ■ b) = gs [ dt(U ■ b) p \ t 

Jo 


t T 

5 

5 S 


-oO(a-I) 


(FI) 
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and 


= gs 2 f d/3 ( da(d fJ 'd IJ ,U)(a(a)) [9 (a — 0) — aO(a — 0)\ 

J Jo 

= gs 2 f dtp [ da(d 2 U)(a(a))(l - a)9(a - 0) 

J o Jo 

= gs 2 f da(d 2 U)(a(a))a(l — a) 

Jo 


N [gU} 2 = Tr(II -b-U-b) 


2 2 
= g s 


dt f dr f da f d/3H fl i y (cr(a))H^ l/ (a(/3)) 
/o Jo Jo Jo 


Ola — max 


o 

t T 

i 

s s 


2 4 

= g s 


-ae(a-j) 

9^0- max - 08 ^0 — - 

J da J dp j d'y J dSH^(a(a))H 111 ' (cr(/3)) 

[0(a — max( 7 , 5)) — a9(a — 5)] [9(0 — max( 7 , 5)) — 09(0 — 7 )] 


= g 2 s 4 f da f d/3 f d'y f d6H^(a(a))H^ v (a(0)) 
Jo Jo Jo Jo 
[min(a, 0) 2 — 2 a/3 min(a, 0) + a 2 0 2 } 

/*1 pOL 

= 2 g 2 s 4 / da(l — a) 2 Tt^(cr(a)) / d00 2 TL^ v (a(0)) 
Jo Jo 

where tools from Appendix ^ were utilized. 


(F2) 


(F3) 


APPENDIX G: Evaluation of N[gU]\ z=0 


In this appendix the ’’loop-generator” N[gU] := — ^Trln(A) for matrix A := 1 — 211 • b for the case z = 0 is 
evaluated. Matrix elements of A are given by 

A^ vt = rf v 8(t - t) - 2(n ■ hy\ T ( = 13 = i |rfj 9 S 2 H^(x)K 0, ^ | (Gl) 

where K(a,0) := 1 — 2max(a,/3) +/3 2 . We will exploit fact that N[gU] can be represented in terms of dimensionless 


variables (a = t/s, etc ...) in the form N[gU] = — ^Trln(A), where A v p = sA ll ’ t=sa ViT=s p = g^ v 5(a — 0) — 
gs 2 'H ,J ’u(x)K(a, 0), where Tr(...) = f* da(...)(a, a). 

The operator K whose kernel is K(a,0) is not completely diagonalisable, since it is not a normal operator, i.e. 

it does not satisfy condition of normality K^K = KK '. To see this, evaluate matrix elements of both sides: 

(K^K)(a, 7 )= f d0K(0, a)K(0, 7 ) = i — a 2 — 7 2 + 27 a 2 + — a 2 7 2 (G2) 

Jo 3 3 

(A'A' t )(a, 7 ) = f d0K(a,0)K(j,0) = i - a 2 - 7 2 + 27 a 2 + - ^a 4 - iy 4 (G3) 

Jo 5 3 6 6 


where 7 > a is assumed, and (|J 14| )-( |J17|) were used. 

We will here hrst resolve the decomposition problem of AT, and only then we will return to evaluation of spectral 
problem for operator A. 


1. Structure of K 

Since K is not normal, it cannot be completely diagonalized. Even worst, its left- and right-hand-side eigenvectors 
are not the same, and those sets of eigenvectors are not complete basis sets. 
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r 

JO 


Here, we will find the right-hand-side-eigenvectors of K. To find them, 
d/3K(a, /3)R(f3) = XR(a) and its two differential consequences in explicit forms: 


start with the eigen-problem 


rl pet 

/ dpRifi) -2a d/3R(/3) - 2 

Jo Jo 


d/3/3R(/3) 


d/3/3 2 R(P) = XR(a) 


(G4) 


- 2 


d/3R(/3) = X R\a) 


(G5) 


- 2 R{a) = XR"{a) (G6) 

Substituting the general solution R(a) = Ae lva + Be - ' 1 '" of that last equation (where X = -^) into ( fd5| ) one gets 
B = A, so R(a) = 2Acos{va). Equation ( |G4| ) further constrains the solution to frequencies v such that _ q 

This gives possible eigen-frequencies v n = mr (n > 1), corresponding to the eigenvalues X n = 2 $ and normalized 

(in terms of scalar product < a\b >:= f* daa(a)b(a)) eigenvectors R„(a) = y/2 cos(nna). Notice that the R n are not 
the left-hand-side-eigenvectors of K. 

Also, although the vectors R n are mutually orthogonal, they are not a complete set: the supposed decomposition 
of unity gives the projection operator Rn{ot)R n ((3) = 8{a — /3) — 1 =: n(a,/3) with defect one. 

A similar procedure for left-hand-side-eigenvectors does not allow solutions in terms of elementary functions, and 
will not be presented here. The only left-hand-side-eigenvector that is easy to obtain is Lja) = 1, with corresponding 
eigenvalue 0. It is not the right-hand-side-eigenvector of K. 

In terms of ortho-normal basis {Lq = 1, R n = y/2 cos(mra)] n > 1} one can write the resolution K = Ko + Kj on 
K D . the diagonalisable part, and Kj , the nondiagonalisable (Jordan type) part: 

OO 

KD = ^ ^ | Rn ^ Rn \ 

n= 1 

oo 

Kj = -^2^(-) n A n |i? n >< Lol =: \KL 0 >< L 0 \ (GT) 

n= 1 


where X n = and \KL 0 >:= K\L 0 >= -V2Y^™=i(~) n ^n\Rn > (note that < L 0 \KL 0 >= 0). 
In the above basis, the matrix form of full operator K is 


K = 


/ 

0 

0 

0 

0 

0 


V 2 X 1 

Ai 

0 

0 

0 


—V 2 X 2 

0 

X 2 

0 

0 



V 2 X 3 

0 

0 

A 3 

0 



-V 2 X 4 

0 

0 

0 

a 4 


V 








(G8) 


In coordinate form, kernels of the K^ j operators are 
r 1 i 

Kj(a, p) = J dyA'(a,7) = - - a 2 

f 1 2 00 4 

Kr>(ct,P)= / d'y K(a, 7)71(7,/?) = - — 2max(a,/3) + a 2 +/3 2 = cos(mra) cos(n7r/3) 

Jo 3 n n ~ 

Also: I\L 0 (a) = | — a 2 . 


(G9) 


n=1 
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2. Spectral structure of A; Evaluation of N[gU\ 2 | 2=0 


A: 


Using the known decomposition of the operator A', one can obtain the corresponding representation of the operator 


A = i (g) I — gs 2 K (£>?{ = 

= |Vo >< L 0 \ + E | R„ >< <g> I - gs 2 > {< R n \ - V2(-) n < L 0 ||j = 

( 1 0 0 0 • • \ 

—gs 2 \f 7 2X\H 1 — gs 2 X{H 0 0 

gs 2 V2X 2 H 0 I - gs 2 X 2 H 0 

~gs 2 V2X 3 H 0 0 l-gs 2 X 3 H 


V 


(G10) 


From that representation one can read its eigenvalues: 


Xa,n (x) :=1- gs 2 A a (x)X n = 1 - 


2gs 2 JX a {x) _ 
7 j 2 n 2 ’ 


n > 1 


where A a (a:) (a = 1 ,D) are the eigenvalues of the matrix H(x). The n = 0 th eigenvalue is equal 1. 
In this way, N[gU] | z=0 is given as 


D oo 


\ T \ rni ( i 2^,s i2 A a (^)\ 1 sm 

= -oEE ln 1 - ~ 2~2 =- ? Tr L ln 


(sy/Wl)\ 


a= 1 n=l x 7 

To obtain the expansion in s, Taylor expand the logarithm in the first form. The first term is: 


(Gil) 


(G12) 


"MiU'-Jee (- 9d ^)=4 fx>«) ^= 4 ^™ 


\a—l 


\n—l 


C i(») 


s 2 


(G13) 


_J_ 7T^_] r _v'-D 

’ A-/71=1 71 2 6 

The second term in that expansion gives 


where S 2 := J2n=i 7? = if and i A a0*0 = Tl 't^(^) = ( d 2 U)(x ). 


=-if: t (4) f-4 ^ 1 ) 2 =- £*■<*> 

a— 1 n=l 


90 


(G14) 


where C 2 {x) := Xf=i A a(z) 2 = Tr L [77(a:) 2 ] = and S4 := J2n =1 F* = fjj- 

Similarly, one can evaluate higher terms in N[gU]\ z=Q . to obtain 


N[gU ]\ z =0 = 1 E ^ 4 (8ffs2)m£m(:c) 

771=1 7 


gs 2 r , s g 2 s 4 4g 3 s 6 g 4 s 8 16g 5 s 10 

- y £tW + -go"A(*) + 2835-A(*) + ^AOr) + + 


(G15) 


where C m (x) := Trz,[W(x) m ]. 


16 






APPENDIX H: On evaluation of \[gU] 


Starting from 


a^(t)=gs daW^aia)) 

J t/s 

b^{ti,t 2 )=gs daH^(a(a)) 

Mt.I) 


and n(fi,t 2 ) = d{t\ — t 2 ) — we obtain: 


(a±)^(t) : = ( a ' n) M (t) =9^ daW^a^a)) 


/' 


(b±±)^(t 1 ,t 2 ) := (n • b • n) Ml/ (ti,t 2 ) =gs daH^u(cr(a)) 


f 


9 [ a -) — a 


ti 


9 I a -- | — a 


9 {a — — | — a 


Then 


9 [ a —- ) — a 


Xo = i [ a±- a±=i f dtg 2 s 2 f da ( dpW ll {a(a))W IJ '(a(p)) 

J J o Jo Jo 

= ig 2 s 3 f da f dp f d"/W ll ,(a(a))W ,l (a(P)) [6 (a - 7 ) - a] [6 (P - 7 ) - p] = 
Jo Jo Jo 

p\ 2 . 

zp 2 s 3 / dot d/3W fX (a(a))W^(a(P)) [min(a,/3) — a/3] = 

Jo Jo 

p 1 /*o: 

= 2ig 2 s 3 da dpW ll (a(a))W ll (a(P)) [min(a, /3) — a/3] 

Jo Jo 

On the same way: 

ii J J a± ■ b±± • ax = 

PS PS P 1 P 1 P 1 

2i / dti / dt 2 g 3 s 3 / da dp dyW, J ,(a(a))Tt' J ' ,/ (a(x))W v (a(p)) 

Jo Jo Jo Jo Jo 


9[P--)-P 


= 2 i 


8 l a — — ] — a 


e[p- t ±)-p 


9 l 7 — max 


ti 


(J19),(718l 


2ig 3 s 5 f da f dp f d'yWi_ l (a(a))H ll ‘'(a(-y))W l ,(a(P))- 
^ Jo Jo Jo 

• / d[i / dz/ [0 (a — //) — a] [0 (/3 — z^) — / 3 ] 0 (7 — max(/z, z/)) = 

Jo Jo 1 

2 ig 3 s 5 f da f dp f dxW ll {a{a))H llv {a(^))W v (a{P))(m:m(a, 7 ) - Q 7 )(min(/ 3 , 7 ) - /Fy) 
Jo Jo Jo 


(HI) 


(H 2 ) 


(H 3 ) 


(H 4 ) 


At c = 0 the a,(3 and 7 integrals can be performed explicitely, using formulas from Appendix [J[ to obtain same 
results as in Appendix|. However, derivation of y m ’s (at 2 = 0) given there is superior to this one (much less tedious). 


APPENDIX I: Evaluation of xVjU)\ z=Q 


From ( G10 ) is easy to obtain the inverse of operator A as: 

OO 

A 1 = \L 0 >< L 0 \ 8 1 + ^2 | Rn >< Rn | <8 (1 — gs 2 X n H)~ 1 + 

n =1 

00 

- V2(-) n gs 2 \ n \R n >< Lq\ 8 (I - gs 2 \ n H)- 1 H 


n —1 


To evaluate y, only projection of that inverse is necessary: 


OO 


A 1 • 7T = ^ \ R ri >< Rn\ 8 (I - gS 2 \ n H) 1 

n =1 


(ii) 


( 12 ) 
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where 7r = X] n =i l-^n >< Rn \ 0 1- To obtain A 1 II multiply above result with s. Then 

p p pi pi °° 

x\ z =o = * / / a ’ A -1 • II • a = is da d(3 R n (a)R n (13)^(sa) 

J J Jo Jo n=1 


1 — gs 2 X n H(x) 


a?(s(i) (13) 


/IV 


where a M (sa) = gs( 1 — a)W //i (a;), W^(a;) := d^U(x). Then 


xl,=o = 2«V E 


n—1 
«.3„2 00 


cio!(l — a) cos(mra) 


L-/o 


W • (I - gs 2 X n H{x))- 1 ■ W = 


7T^ ^ TIT* 


(14) 


From this point one can proceed in two alternate routs. First is to Taylor-expand matrix factor in the above form 
and to obtain an explicit expressions for \m s (x = Y^m-o Xm)' 


where 


Xm\ z —o — 


2m+1^2m+3^m+2 
^■2m+4 


(W ■ H m ■ W)S 2 m +4 


° n / 'inl2 9 /" 92 m +4 _ i \ 2m+4 

■W* = E JJJ = 4 (1 - 2- 2m - 4 ) C(2 m + 4) = ( -|B 2m+4 | 


n=l 


(2m + 4)! 


(15) 


(16) 


giving 


2 m+2 (2 2m + 4 - 1) s 2m+3 3 m+2 |B 2m+ 4| 

XmL=o = *---- 1 2 +4 {W ■ H m ■ W) 


Several first members of that series are: 

Xol—o 

Xi I,=o = 


(2m + 4)! 

i s 3 g 2 W 2 
12 

i s 5 g 3 W -H-W 


XsU =0 = 


The second possible route is to regroup terms in 


L~=° 6Q 

i 17 s 7 g 4 W ■ H 2 ■ W 

X2 ' z =° ~ 5040 

i 31 s 9 g 5 W-n 3 -W 


37800 




71=1 


where a 2 := ^§- H , and to evaluate given sum 

' [1-H n ] 2 1 


Then 


xl z= o = 


E 

71=1 


|£V w 1 f tan(^) 

^ 2 « 2 V (¥) 


7 r 2 / tan (£2) 


2« 2 l (4?) 




-1 


tan s 


\ 


- 1 


W 


(17) 


(18) 


(19) 


( 110 ) 


(Ill) 
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APPENDIX J: Tools for parametrix calculations 


(z.d x + n)k n (x,y) = jn(x,y) ■ 

d?f(a(a)) 

[ data™f(a{a/3)) 
Jo 

f d/3/3 l [ daa m f(a(aj3)) 

Jo Jo 

f 1 d(3(3 n f 1 daf(a(3) 
Jo Jo 

f dyy 1 [ daoT f d/3/3 n f(ay)g(/3y) 

Jo Jo Jo 

f dyy 1 [ daa m f d/3(3 n f(ay)g((3y) 

Jo Jo Jo 

[ daf(a) ( d(3g({3) 

JO J a 

f d/30(a — max(/3, 7 )) 
Jo 

(n ■ b)^(ti,t2) 
(n • 6) MI /(ti, t 2 )| 2 _ 0 

[ dy f dS9(a — max( 7 , 5))9{/3 — max( 7 , S)) 
Jo Jo 

[ dy f d89(a — max( 7 , 5))9((3 — 7 ) 
Jo Jo 

[ d/3/3 n max(/3, a) 

Jo 

/ d/3 max(/3, a 
Jo 

[ d(3/3 2 max(/3, a) 

Jo 

/ d/3 max(/3, a) max(/3, 7 ) 
Jo 

f dy [9 (a - 7 ) - a] [9 (J3 - 7 ) - /3] 
Jo 

9 (7 — max(/i, v )) 

f d/3(min(a,/3) — a/3) 
Jo 


kn(x,y)= I daa n 1 j n (a(a),y) 


L 

a m (df)(a(a)) 


p-™- 1 ( daa" 1 f(a(a)) 
Jo 


l- 

1 


~m L doi ^ am ~ 

f da(l — a n )/(a) 
Jo 


1 


l — m — n — 1 


f 


daa™ /(a) < (1 — a 


+ f 1 d/3/3 m (l — f3 l ~ m ~ n ~ 1 )g(f3) 
J a 


) [ d/3/3 n <?(/3) 
Jo 


l — m — n — 


gs 


= gs 


- -r [ da(a m -a l n 1 )f(a) [ d(3/3 n g((3) 

-n-1 J 0 Jo 

[ dag(a) [ d/3f((3) 

Jo Jo 

a 0 (a — 7 ) 

f da (d^dv U) (a (a)) — f daa(d IJ .d v U)(a(a)) 

max(ti,t2)/s Jt^/s 


1 / tl t2 \ 1 ( t2 

__ max 


H^{x) = i gsTl flv (x)K ^ fl tz 


= min(a, /3 ) 2 

= amin(a,/3) 

1 


(n + 2) 


1 + 


n+2 


n + 1 


1 


) = ^(l+«) 




11 . . , 1 . .3 

- + -7a mm(7, a) + - max(7, a) 


— 7) — /3] = min(a, /3) — a/3 

= # (7 - M) 0 (7 - ") 
= |a(l-a) 


(Jl) 

(J2) 

(J3) 

(J4) 

(J5) 

(J6) 
(J7) 
(J8) 
(J9) 
(J10) 
(Jll) 
(J12) 
(J13) 
(J14) 
(J15) 
(J16) 
(J17) 

(J18) 

(J19) 

(J20) 


APPENDIX K: Schwinger-DeWitt calculations 


Rewrite the system (9.4) in the form 


{z- d z + n)k n [y, z ] = j n [y, z] 


(Kl) 


where n > 1, k n [y,z] := k n (x,y) and: j i = gU, j 2 = ~d z k 1: j 3 = -d 2 z k 2 - (3+i) 2 , j 4 = -d 2 k 3 - 2{d z k^\{d z k,2}, ■ ■ ■■ 
Then, for every given (known) j n \y, z], the corresponding k n [y, z] is given by the expression (Appendix [j], Eq([Il|)) 
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kn[y,z\ 


(K2) 


daa n 1 j n [y,az] 

Then, with known k n [y,z], one can evaluate the next source function j n +i[y,z], etc ... During the calculations 
several characteristic procedures are taking place. As an illustration, let us calculate the first several ks. 

Since ji[y, z] = gU(y + 2 ), it follows that ki[y, z] = g f 0 daU(y + az). 

The second source function takes the form j 2 [y, z\ = —gd 2 f ( ' daUfy + az). Then, each d z converts into ad y upon 
action on U(y + az). So j 2 [y, z] = —g fg daa 2 (d 2 U)(y + az). That was first characteristic procedure: d z —> ad y . 
Then k 2 [y, z\ = fg daaj 2 [y , az] = — g fg daa fg df3(3 2 (d 2 U)(y + (3az). Rescale (3 — »• /3/a to make argument of U free 

of a: k 2 [y, z] = —g fg daa fg d(3(3 2 a~ 3 (d 2 U)(y + (3z) = —g fg daa~ 2 ff d/3/3 2 (d 2 U)(y + f3z) (The second procedure: 
/3 —> /3/a). Extend the range of integration of (3 back to [0,1], introducing the step function 9(a — (3 ), switch the 
order of integration (first to perform the a integral, and only then the [3 integral), and limit the range of integration 
of a to [/3,1] (according to the argument of the 9 function): k 2 [y,z] = —g fg d(3f3 2 (d 2 U){y + (3z) f^daa~ 2 (Third 
procedure: switch the order of integrals). The integral over a can be performed, giving daa~ 2 = /3 -1 (l — (3) and 
fa[y,z] = -gfg d(3(3{l - f3)(d 2 U){y + (3z). 

The third source function becomes js[y, z] = —d 2 k 2 — {d z k\) 2 = g fg d/3f3 3 (l ~ f3)(d 4 U)(y + f3z) — g 2 fg d/3(3(dU)(y + 
(3z) fg dyy(dU)(y+'Yz). Replaying the same order of procedures as above, one can obtain (in seven lines of derivation) 

* 3 [y, z\ = f fg daa 2 {l - a) 2 (d 4 U)(a(a)) - 2 g 2 fg da( 1 - a)(dU)(a(a)) /“ df3(3{dU){a((3)). 

In the same way, the fourth k function can be obtained as: 

1 r 1 

ki[y,z] = —g J daa 3 ( 1 - a) 2 (d 6 U)(a(a))+ 

pi pOi 

+g 2 daa(l — a) 2 (d 3 U)(a(a)) / d/3/3(dU)(cr{/3))+ 

Jo Jo 

pi pcx. 

+2 g 2 / da(l - a) 2 (d fl d l ,U)(a(a)) / d/3/3 2 (d^d^U)(a(l3))-\- 

J 0 Jo 

+g 2 [ — (1 - a) 2 (dU)(a(a)) f d(3f3 3 (d 3 U)(a(f3))+ 

+g 2 [ —(1 -a)(dU)(cr(a)) [ d/3/3 2 [2a - /3(1 + a)\(d 3 U)(cr{f3))+ 

J 0 a Jo 

+g 2 [ daa{dU)(a(a)) [ d/3/3(l - /3) 2 (d 3 U)(a{(3)) 

Jo Ja 

(K3) 

This takes about 30 lines (one and one-half pages of calculation), that is a factor of four longer than the evaluation 
of k 3 . One can expect that the length of computation for higher order functions will grow at least as a geometric 
progression. It is interesting that all steps (procedures) are straightforward, and one could be tempted to program 
some symbolic mathematical tool to perform them. 

Some of the time-saving formulas obtained during calculations are placed in Appendix |jj. 

It is interesting that if one asks different kind of questions, then certain exact (i.e. in all orders in s) partial 
information can be obtained. For example, in the Appendix |] the first order in g contribution to k is derived. 



APPENDIX L: Derivation of Fi(z|s) 

From (|9.2|) follows the IVP for Fi(z|s): 

[-i(sc?g + z ■ d z ) + sdl\ Fi(^|s) = sU(x) 

F!(z\0) = 0 (LI) 

Then, using Schwinger’s parametrization one can write 

pOO 

F 1 (z\s) = i dte-^-^+^+^lsUix) (L2) 

Jo 
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For a moment switch to an algebraic structure by introducing differential generators a := —sd s , b := —z ■ d z and 
c := —isd 2 with the commutational relations [a, b] = 0, [a, c] = — c, [a, b] = 2c. Then we are interested in r := e t ( a + h + c ). 
Assume that it can be factorized in form f = e A (t)a e B(t)b e C(t)c_ Then the (|lT[) gives equation 


a + b + c = Aa + Be A Ad (a) 6 + Ce A Ad (a) e B Ad (S) c 

Since e A Ad ^b = b, e B Ad = e 2B c and e A Ad ^c = e~ A c, the IVPs for A, B and C follow: A = B = 1, 
C = e A ~ 2B , A(0) = 5(0) = C(0) = 0. Their solution is A(t) = B{t) = t, C{t) = l-e~ t giving: f = e ta e tS e(r- e_ ‘) e . 


So: 


F 1 (z\s) = i / dte- tsd °e- tzd *e- is ^- e ^ d2 *sU(x) 


Use now the Fourier representation for U{x ) {x = y + z) to obtain: 

F 1 (z\s) = j dpe ipv i dte- tsds e- tz ^e ipz e +ia ^- e ~ t)p \u{p) 

= J dpe ipv i dte-^’e^-* e +is ( 1 - e ~ t)p2 sU(p) 

= J dpe ipv i J dte^'e +ise ~ t{1 - e ~ t)p2 e^sUip) 

= isfdte-* j dpe i A v+ze ~ t )+^~ t ^~ t )p 2 U{ p ) 


i.e. 


pOO 00 -j 

FUzIs) = is / dte~ t 'S~' —- (—*se _t (l — e~ t )d 2 ) n U (y + ze - *) 

Jo ^o nl 


Introduce a = e 1 and the final expression is established: 

(*1 OO 


Z* 1 ^ 1 

FUUs) = is / da — (—isa(l — a)d 2 ) n U(y + za) 

/n n! u 

1/0 n=o 

= is f daU(a(a))+s 2 f daa(l — a)d 2 U(a{a))+ 

Jo Jo 

— — [ daa 2 (l — a) 2 d 4 U(a(a)) — — [ daa 3 {l — a) 3 d 6 U(a(a)) + ■ 
2 Jo 6 Jo 


(L3) 


(L4) 


(L5) 


(L6) 


It is easy to check that first four terms indeed do agree with expressions obtained through more tedious calculations 
in previous Appendix. 

At ^ = 0 the a dependence disappears from argument of U and integrals over a can be performed (f Q daa n (l—a) n = 
B(n + l,n + l) = (2 " ! "- )! ), giving 


Zi(~|0) = is^T 


n —0 


(2n + 1)!! V 2 


-a2) U(y) 


= isU{y) + j(d 2 U)(y) l -^(d 4 U)(y) - ^_( d e U)(y ) + 


(L7) 


Julian Schwinger, On Gauge Invariance and Vacuum Polarization; Physical Review, Vol. 82, P.664-679 (1951) 
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